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ABSTRACT 

The 'Veathering vaning" motion of a wind turbine with a moving rotor is an oscillatory problem with 
a variable moment of inertia. The analysis of such a motion rec^uires the solution of a non-linear 
differential equation. In this article an approximation method is presented for reducing the problem to 
an equivalent constant moment of inertia problem. 

The method is based on the assumption that a moving rotor is an integrator and, therefore, the 
problem will behave as if it has an averaged moment of inertia. It is further assumed that this will be 
a valid solution to the problem if the rotating speed of the wind turbine is infinite. The method 
consists of determing the integrated average of the moment of inertia for a single rotation. This 
averaged value can then be used to determine equivalent natural frequency of the sy¥fem and other dynamic 
propert ie s . 

The method is shown to be valid by solving the non-linear differential equation for various rotating 
speeds. It was found that the cycle time is the equivalent cycle time if the rotating speed is 4 times 
greater than the systems minimum natural frequency. The ratio of equivalent to minimum cycle time is 


t/to - Clmax/Imin)^ . 


INTRODUCTION 

The oscillatory motion of a system with a 
variable moment of inertia 1$ complex. The 
analysis of such a motion requires the solution of 
a non- linear differential equation. In this 
article a method for reducing the problem to an 
equivalent constant moment of inertia problem is 
pro rented • 

The application which generated interest in 
this problem is the "weather vaning” of a wind 
turbine with a moving rotor. The method, estab- 
li sited here is applicable to any problem in which 
the moment of inertia is a variable about the 
axis of rotation such the roll of a helicopter. 

The problem is illustrated in Fig. 1, It is 
desired to describe the motion of the system 
about the axis A- A while the rotor is moving. One 
caiT^^ from Fig. 1 that the moment of inertia 
about A-A varies from a minimum when the rotor 
is vertical (Position A) to a maximum when the 
rotor is horizontal (Position B). This variation 
is continuous and cyclic with each turn of the 
rotor . 

One can also see from Fig. 1 that the problem 
cannot be handled with a single initial condition. 
A different and unique motion is obtained for each 
initial position of the rotor. Hence, the problem 
has a stochastic nature. 

= MALYSIS 

At any rotor position, ibe instantaneous 
moment of inertia of the rotor, I*, about axis A-A 
is : 

2 2 2 
I. — mr, sin i = I sin di 
X k ^ r ^ 

where m is the mass of the rotor. The terms r^^ 



Fig. 1 - Problem Geometry 


and are the radius of gyration and the moment 

of inertia of the rotor about its axis of rotation 
respectively. The moment of inertia of the rest 
of the system about A-A, which is constant, is: 

I = I + mh^ 
o c 

where I is the moment of inertia of the center- 
body anS h is the distance from the center of the 
rotor to axis A-A. The expression for the total 
moment of inertia I is the sum of the constant: and 
variable portions 

I = I + I. = I + I sin^^i 
o 1 o r ^ 
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or 


I - ^qCI Jsin (1) 

where J = I /I . 

r o 

When the entire system is oscillating about 
axis A-A with a moving rotor, it is apparent that 
the rotor is acting as an integrator of moment of 
inertia. It is hypothesized that if the rotor has 
an infinite angular velocity the moment of 
inertia I behaves as a constant integrated aver- 
age value. It is this hypothesis which is the 
basis of method presented here. 

To establish the utility of the method, 
assuming for the moment that the hypothesis is 
correct, there are two questions which must be 
answered : 


1} What averaging technique should be used? 

2) How close must the speed be to infinite 
for the method to be useful? In other 
words, how fast is fast? 

The hypothesis is proven by examining the 
simple harmonic motion of the system in Fig, 1. 
The proof is valid for more complex situations 
such as those including damping and forcing 
functions. In demonstrating the proof, the 
method of solving a variable moment of inertia 
problem will be established. 


The equation for simple harmonic motion is: 


9 2 
d g _ -Kg _ ^ ^0 ^ 

dt^ I 1+Jsin^(wt) 


( 2 ) 


and 


= K/Io 


(3) 


where K is the torsional constant and w is 

o 

defined in Eo. 3 as the reference freauencv for 
the system. In writing Eq. 2. the moment of 
inertia I has been replaced by Eq. 1 and the 
rotor angle 4) has been replaced by rotor angular 
velocity w and time t. 


Eq. 2 is dimensional ized bv letting 


T t/t 

0 

n = oj/w 


w t 
o 


and 


0 = a/g 


where g^ is the initial angular displacement. 

Substituting these expressions into Eq. 2, one 
has 


d^e 

dx2 


2 

1+Jsin (nx) 


- -ke 
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The term k is defined as the natural frequen- 
cy ratio 

- (0 2 

k = i2 = ) C5) 

1+Jsin (nx) w 

0 


and can be used to obtain the natural frequency of 
the system at anytime. 


The hypothesis states that as the rotor angu- 
lar velocity approaches infinity, the natural 
frequency ratio k approaches a constant integrated 
average value. The average value is obtained by 
integrating Eq. 5 over a quarter cycle 


k 




d4> 

2 

1+Jsin 4) 


7t/2 

2 r d4> 

7T V .T - 2 2 

o (J+lJsin 4>+cos 


C6) 


From a table of definite integrals*, one obtains 


7T/2 

f 


o 


dx 


2 2 2 
a sin x+b cos 


2 

X 


TT 

2ab 


k 




(7) 


Eq. 7 is the answer to the first question. Since 

J + 1 is I „^/I 4 ^, a more convenient form for 
max min 

Eq. 7 is 


t/t = Cl /I ^ ) 
o max min 


( 8 ) 


The hypothesis is proven if it can be shown that 
the system oscillates at a cycle time defined by 
Eq, 7 as n increases. It should be noted in 
passing that the averaging method defined by Eq. 6 
is the proper one. Initially, the following 
averaging was used which is incorrect. 




1+Jsin (|))d(|i 


The error the author made was that reciprocal of 
the average is not the average of the reciprocals. 

Eq. 4 is a non-linear differential equation. 
It was solved numerically for a number of cases of 
n and J using the Continuous System Modeling 
Program (CSMP) which is standard IBM software. 

The solution to thp differential equation 
with time ratio x is presented in Fig. 2 for J = 1 
One can see that the frequency shifted signifi- 
cantly with increase in rotor speed from n = 0 to 
n = 1. Fig. 2 also shows a significant change in 
displacement with speed ratio n. In 3 reference 
frequency cycles, the solutions for other than n = 

0 does not appear to be repeating indicating the 
stochastic nature of the problem. 


*Handbook of Chemistry and Physics, Chemical 
Rubber Publishing Co. 
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Fig, 2 - Displacement Response with 
Variable Rotor Speed 


In Fig. 3 the time ratio t for the first 
cycle is presented with increasing rotor speed 
ratio h7~” The limiting value as pfed^icted by Eq, 
7 is Indicated in Fig. 3. One can see from Fig. 
3 that the time ratio equals the limiting value 
when n > 4. Hence, the hypothesis is proven and 
the second question is answered. 


The curves in Fig, 3 are smooth even though 
the Fig. 3 does not indicate it. The curves all 
start at T = 2tt because the problem was always 
started with the system at minimum moment of 
inertia or with rotor at Position A. If the prob- 
lem were started at maximum moiuont of inertia, the 
solution would have started at ^ 

At any other position the solution would start 
anywhere in between; It Tsthe variation in 
moment of inertia with rotor position that causes 
the oscillation of the first cycle time ratio in 
Fig. 3 before It "damps” down to the limiting 
value , 

It appears from Fig. 3 that as J increases 
the first cycle time ratio approaches the limiting 
value more quickly. The reverse is also true. 

The solution for J = 0 is a constant horizontal 
line as indicated in Fig. 3. The speed ratio for 
utilizing the limit solution may maximijinize with 
J. 


CONCLUDING REMARKS 


In conclusion, it can be stated that: 

1) The hypothesis is true. 

2} A variable moment of inertia problem can 
be analyzed as a constant moment of inertia prob- 
lem if the rotor speed is 4 times the reference 
frequency, which is based on minimum moment of 
inertia. 

3) The method of solution is to multiply 
the reference frequency by the fourth root of the 
minimum to maximum moment of inertia ratio to 
obtain the system’s natural frequency. 

4} If n > 4 the initial starting point is 
not important and the problem is deterministic. 
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Fig. 3 - First Cycle Time Response 



QUESTIONS AND ANSWERS 


E.W. Beans 


From: T. Base 

Q: When you set up your initial equations, why didn't you write them in form 

^ (I<^) = Applied Torques {including damping) 
so that: ^ = Applied Torques 

then the rate of change of with time could be used directly in the equations. 

A; The equation you wrote is not the one I wanted to solve. Your equation appears to 

have a velocity dependent term. Since dl/dt varies cyclic I would try to solve it 
using the averaging technique , 

From: Art Smith 

Q; Do you plan to check your results with a more exact method? 

A: 1^0. My results appear to satisfy application for which it was developed. 

From: G. Beaulieu 

Q: You have considered rigid blade for your analysis; could the harmonic deformation 

of a vibrating blade significantly change the moment of inertia? 

A: YeSj if the mass center is displaced. 

From; Dan Schiff 

Q: Does your solution account for dynamic effects, e.g., gyroscopic effects — or only 

static? 

A; ffo, but the method should he applicable to the solution of any equation with cyclic 
coefficients , 
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